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1. Introduction 

Modern calculation physics has several important features, which may push ahead this part of 
science. To our point of view among them are: 

1. the requirement of reliability of the complicated numerical or analytical calculations, related 
to theelimination of bugs in the computer codes and final results. It seems to be obvious to 
realize this requirement, but sometimes not simple (see e.g. [[I]]); 

2. the requirement to avoid shortcomings of definite numerical calculations, which may mani- 
fest themselves in the non-proper rounding off errors in the numerical results (see e.g. [JT|] , 

0D; 

3. the interest to phenomenological applications of various perturbative results and the devel- 
opment of different resummations approaches of perturbative series (see e.g. [0], [Q]); 

4. the appearance of strong intention to go beyond straightforward theoretical calculations and 
the desire to find new mathematical [||] and physical [|6|] properties, related to the differences 
in the topological structure of various multiloop Feynman diagram, which define special 
analytical structures of perturbative series under investigation. 

I will try to demonstrate that continuing at present QCD calculations of characteristics of deep- 
inelastic scattering (DIS) satisfy all requirements mentioned above. In particular, I will concen- 
trate myself on discussions of possible relations of the concrete analytical structures in the con- 
crete results of perturbative QCD calculations to the existence in nature of definitely speaking 
non-perturbative effects, namely light quarks flavour asymmetry u(x) < d(x) (see Ref.[^]) Related 
phenomenological manifestations of this effects will be also considered. Presented in Ref . [||] con- 
jecture about existence of differences between QCD corrections to the non-singlet (NS) moments 
of charged leptons - nucleon and neutrino-nucleon DIS, which manifest themselves in the large 
A^-limit, proposed in Ref.[||], will be described.The perturbative indications to the generation of 
s(x) — s(x) strange quark asymmetry (see Ref.[{7]] will be mentioned as well. 



2. Definitions of some perturbative series in DIS 



Let me follow normalization conditions of the related perturbative series, used in the work of 
Ref. [Q] devoted to the detailed fits of the experimental data for the NS xFi, SF of vN DIS, provided 
by CCFR collaboration [|l(i|]. Consider DGLAP equation in the NS case 



Q- 



dQ< 



dz 



PNs(z,a s )+P(a s 



dlnC NS (z,a s ) 
da s 



7 NS 



X 

z 



(2.1) 



where Cns(z,ci s ) and Pns(z,ci s ) are the NS combinations of the coefficient functions and splitting 



functions, F 1Vit - NS combinations of F 2 / v structure functions (SFs) of charged leptons ( I) or V - 
nucleon DIS. The QCD j8 -function is defined as 



n>0 



n+2 



(2.2) 
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where a s = a s / (Ak) and the first two scheme-independent coefficients are normalized as 

A) = (y Ca - |jV>J (2.3) 
34 , 10 \ 

—C\- 2C F N F - —C A N F J (2.4) 

The general expressions for the Casimir operators of the SU (N c ) group of colour are Ca = N c , 
Cf = (Nc — l)/(2N c ), while Np is the number of quarks flavours. The series for the coefficient 
functions and splitting functions of Eq.(2T) read 

C NS (z,a s ) = £cf(z)< 



n>0 



P N s(z,a s ) = £i&K 



ji+i 



(2.5) 
(2.6) 



n>0 



The corresponding non-expanded in a s NS splitting functions we will be interested in have the 



form, given in the work of Ref . [11] 



mk 



P-- — S rP V + P S 
Pqtqk = SikPqq + Pqq 



(2.8) 
(2.9) 



where the next-to-next-to-leading order (NNLO) o? s corrections to the NS combinations of valence 
( Pq q , P^) and sea (P^(jc), P^) contributions were derived analytically in [11]. 

These results are related to the coefficients of the anomalous dimensions of NS moments 
through the following Mellin transformation 



+ 1 



n>0 



dx 



P NS (x,a s ) 



^- l = -2f l dxWp^{x) 

JQ n>0 



x*- 1 (2.10) 



In the case of N=2,4,6,8 Mellin NS moments of the charged lepton-nucleon DIS SF F2 they were 
previously calculated analytically at the NNLO in Ref. fll8|], N=10 NNLO results were obtained in 
Ref.Jl9|], while for N=12,14 they are known from the calculations of Ref.[pO|]. In the case of N=16 
NS moment of F2 SF the direct calculation of the 0(a%) approximation for yS s (a s ), performed in 
Ref. [gj]|, helped to check the three-loop results for P^s(x,a s ) [|TJ] by verifying coincidence of 
anomalous dimension of N=16 moment, calculated in Ref.[j2~H], with the result obtained from the 
expression of Ref. [11] using Mellin transformation in Eq.( |2. 10 ). 

It should be stressed, that among the first phenomenological applications of the results of 
Refs.[0],[|l|] were the NNLO DIS fits of Refs.[g2|],[|l^], [0],[||],@, @, where the works of 
Refs,[||], [f|] were devoted to the analysis of F 2 SF data. The works of Refs.[||],[0], [|e]],[[T7|] 



were aimed to the "approximate" NNLO fits of CCFRXF3 data for V./V DIS, which were essentially 
based on the results of rather complicated and, to my point of view, distinguished calculations of 
the NNLO corrections C^ 3 (jc) to the coefficient function of the DGLAP equation for xF 3 SF of V./V 



DIS, performed by van Neerven and Zijlstra [22]. The first confirmation of this result came from 
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the coincidence of its first Mellin moment with the NLO analytical approximation for the Gross- 
Llewellyn Smith sum rule, evaluated by Gorishny and Larin [23]. More detailed cross-checks were 
done in Ref.[|4|] using different methods. 

The concrete theoretical arguments, that for moderate and high N the NNLO approximations 
of the NS anomalous dimensions for F2 SF will be rather closed to the NNLO approximations of NS 
anomalous dimensions for xF-$ SF were used in Ref. [13]. Theoretical studies of these arguments, 
which in part are based on the proposal to use extrapolation procedure from even to odd values of 
N [ 12], are on the agenda (careful checks of the validity of the extrapolations procedures was done 
in Ref.[||] and in the resent work of Ref.[25]). 

The "approximate" feature of the NNLO fits of xFt, data pushed a-head analytical calculations 
of the NNLO approximations for NS anomalous dimensions and coefficient functions of xF$ SF of 
odd Mellin moments with N=l,3,5,7,9,ll,13 @ (the N = 1 NNLO results of Ref.@ coincide 
with the NNLO calculations for the Gross-Llewellyn Smith sum rule [^6|]). 

Let us now return to the definition of the coefficient function in Eq^l^) in the case of NS 
combinations of F2 SF. The NNLO term of order a 2 s can be divided in two parts 



C 



.NS,F 2 



(x)=C^ + \x)±C^-\x) 



(2.11) 



where the sign plus in between two terms corresponds to the NS contribution to the coefficient 
function of the F2 SF of charged lepton-nucleon DIS, while the sign minus appears, when the 
coefficient function of the NS contributions to F2 SF of neutrino-nucleon DIS is considered. The 
expressions for the terms C^'^ix) and C^^-\x), were calculated in Ref.[E7f] and confirmed 
in Ref. [24]. Thus, all discussions presented above, are demonstrating that the first requirement, 
mentioned in the Introduction, i.e. the reliability of the results obtained, is fulfilled in the case of 
the complicated NNLO calculations of DIS characteristics. The third requirement, and in particu- 
lar, its first part, namely the applications of new perturbative results for the phenomenological 
studies, is also fulfilled in the case of DIS characteristics. 

Indeed, the appearance of the calculations of Ref. [20] gave the possibility to perform more 
careful analysis of CCFR xFi, structure functions data, which came from the studies of vN DIS 
process at Tevatron. These phenomenological fits were made in Refs.[29], [0],[3C],[31], [32]) us- 
ing different theoretical approaches. Note, however, that after these works the number of questions 
still remain for possible future studies. For example, in spite of the fact, that the NNLO theoretical 
inputs of various fitting codes of analyzing DIS experimental data ( i.e. the ones, based on applica- 
tion of the DGLAP equation and on the reconstruction of structure functions from their moments 
using Bernstein and Jacobi polynomials) are coordinated, the outputs sometimes differ drastically. 

Among the concrete questions is the existence of "small" NNLO value a s {M z ) = 0, 1 142 ± 
0.0008 obtained by Bernstein polynomial analysis of the CCFR xF^ measurements [|3l|], which is 
agreement with "small" NNLO result a s (M z ) = 0.1141 ±0,0014 of the combined NNLO fits of 
various DIS data [33] (other "small" NNLO values were recently discussed in Ref.[34]), but is 
not consistent with the "world average" result GC S (M Z ) = 0.1189 ±0.0010 [|5|]. It is also worth 
to mention, that, that the "small" result of Ref.[31] does not also agree with other NNLO ap- 
plications of the Bernstein polynomial technique for studies of the same data CCFR xF^ data, 
which give the following results a s (M z ) = 0.1 166 ±0.0013 @, a s {M z ) = 0.1196^;^ and 
0Jj(Mz) = 0.1189 ± 0.0019 [[O]. It will be of interest to trace the origin of these disagreements of 
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applications of the same data and same theoretical method in the different works on the subject. 
Moreover, "small" NNLO a s {Mz) values, extracted from various characteristics of DIS, do not 
agree with the the recent NNLO result a s (M z ) = 0.1 19 ±0.002 (exp) ±0.003 (theory), obtained 
within the updated set of NNLO partem distributions of UK group [06|]. It should be stressed, that 
this value is in perfect agreement with the one, obtained from the NNLO fits of xFj CCFR data 
using Jacobi polynomial technique, namely a s (Mz) = 0.119 ±0.005 (exp) ± 0.003 (theory) [Q], 
where experimental error is dominated by systematic uncertainties of CCFR vN scattering DIS 
data. We think that more careful coordination of theoretical and experimental uncertainties of dif- 
ferent NNLO determinations of a s (Mz) from the characteristics of DIS data will help to resolve 
existing disagreements. 



3. Adler and Gottfried sum rules: the definitions 



Two other sections are following the original presentation at ACAT07 Workshop of the ma- 
terials of the work, made in collaboration with Broadhurst and Maxwell (see Ref.[[J ). Some 
additional considerations [37], which are supporting the results of Ref.[Q] will be also commented. 

The isospin Adler sum rule [38] is expressed through the structure function F% N of neutrino- 
nucleon DIS as 



l dx 
o x 



F^ p (x,Q 2 )-F 2 v "(x.Q 2 



4h 



In terms of partem distributions Eq. (3.1) takes the following form 

Ia = 2 dx[u v (x) — d v (x)] = 2 
Jo 



(3.1) 



(3.2) 



where u v (x) = u(x) — u(x) and d v (x) = d(x) — d(x) are the valence parton distributions. It was 
shown, that Ia receives neither perturbative nor non-perturbative (l/Q 2 ) -corrections [39]. In view 



of this the Adler sum rule is Q independent and demonstrates the consequence of the property 



of scaling [4C]. This property (or so called "automodelling" behavior of structure functions) was 



rigorously proved in the case of charged lepton-nucleon DIS [ |41| ] with application of general prin- 
ciples of local quantum field theory, described e.g. in the classical text-book jp2[]. It should be 
stressed, that the validity of this sum rule is supported by the existing neutrino-nucleon DIS data, 



which show no significant Q variation in the range 2 GeV < Q < 30 GeV and give [ |43| ] 



.exp 
A 



2.02 ±0.40 . 



(3.3) 



Though the error-bars are quite large, the precision could in principle may be improved by analyz- 
ing data of other viV DIS experiments in more detail. 

However, it is known, that in QCD scaling is violated. The sources of its violation are related 
to the asymptotic freedom effects and non-perturbative contributions. Both types of these effects 
manifest themselves (though is some puzzling way) in the analog of the Adler sum rule, namely in 
the Gottfried sum rule. It can be defined as the first N = I NS Mellin moment of the difference of 
F2 SFs of DIS of charged leptons on proton and neutron, namely 

f 1 dx 
h x 



F' p (x,Q 2 ) 



F{ n (x.Q 2 



ijT dx(u v (x,Q 2 )-d v (x,Q 2 )^j = ~ (3.4) 
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The definition of Eq. ( pT4| ) is presented in the case of assumption accepted previously that the sea 
quarks are flavour-independent. It corresponds to the condition u[x, Q 2 ) = d(x, Q 2 ), accepted in the 
early works on the subjects. However, due to the appearance of experimental data for the muon- 
nucleon DIS, Drell-Yan process and semi-inclusive DIS we know at present that this condition is 
violated and u(x, Q 2 ) < d(x, Q 2 ) ( for reviews see, e.g. [44]-[|46|]). Therefore, the definition of the 
Gottfried sum rule was modified as: 



dx 



o x 



F l 2 p {x,Q 2 )-F> n {x.Q 2 



f G + - dx(u(x,Q 2 )-d(x,Q 2 ) 



(3.5) 



where the last term is related to the manifestation of isospin-breaking effects in the Dirac sea. 



Without it the valence part of the Gottfried sum rule of Eq. ( |3.4[ ) contradicts the existing result of 
the analysis of muon-nucleon DIS data, performed by the NMC collaboration, namely [ |47| ] : 



I G (Q 2 = 4 GeV 2 ) = 0.235 ± 0.026 . 



(3.6) 



Note, however, that the error of this experimental result may be over two times larger due to more 



careful treatment of the uncertainty from the small x region [48]. In spite of this, definite contradic- 



tions of theory with the quark-parton result still exist [ }47p and it is worth to study the contributions 
of the perturbative and non-perturbative corrections to the valence part P G . 



4. Large iV c -expansion of the Gottfried sum rule and its relation to the Adler sum 
rule. 

Let us now discuss the results of calculations of perturbative QCD corrections to the valence 
part of the Gottfried sum rule, following the considerations of Ref. ^ and using the ideology 
of large-A^ c limit, proposed in Ref.[j8]]. The solution of the renormalization group equation for the 
valence contribution P G to the Gottfried sum rule has the following form 



■ A(a s )C {l) (a s ) 



with the anomalous-dimension term 

A(a s ) = exp 
1 



(x) 



S P(x 

1 (7i C7V=1) ^ 2 " {N=]] 
2~ 



dx 



1 + 



i ri 



IN=1) 



+ ■ 



A) p 2 A> 

The anomalous dimension function of Af = 1 moment is defined as 



ft 

a s 

K 



(4.1) 



(4.2) 



(I) N=l 

Yns 



(4.3) 



«>0 



The first coefficient yjy N 1 is identically equal to zero [[49TJ . The leading correction to Eq.(^3j) 
comes from the scheme-independent two-loop contribution [50] 



7i 



(0 jv=i 



-4(C 2 -C A C F /2)[13 + 8C(3) - 12C(2)] 



(4.4) 
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confirmed later on in Ref.[]5~ij], Notice the appearance of the distinctive non-planar colour factor 

(Cj-C A C F /2) = 0(N°), which exhibits 0(l/N%) suppression at large-/V c , in comparison with the 

individual weights of planar two-loop diagrams, namely Cp and CfCa, that are canceling in the 

f „,(') Ar=i 
expression for /j ' 

The three-loop coefficients of anomalous dimensions of N-th moments moments for the NS 
combinations of structure functions F 2 of neutrino - nucleon and charged lepton-nucleon DIS are 



related to the Mellin moments from the splitting functions as [11] 



72 



(v)N 



7 2 



2 dx 
Jo 

2 dx 
Jo 



,(2)-, 
NS I 



,(2)H 
- NS 

p(2)H 



x^ 1 
x^ 1 



(4.5) 
(4.6) 



Taking into account the nullification of anomalous dimension function for the Adler sum rule, one 
can derive the useful identity 

^1 



7 2 



.(') N=l 



dxP, 



(2)- 







NS 



(4.7) 



Integrating the results Ref . [|T1|] and taking into account mathematical analytical procedures of 



Ref.[|52[], in Ref.[|6j] the analytical expression for N 1 was obtained: 



7 : "" =1 = (C 2 F -C A C F /2)\c F 



SI) N=\ 



290 - 248^(2) + 656^(3) - 1488^(4) + 832£(5) + 192^(2)^(3) 



+ C A 



1081 + 980 



9 3 
304 

161.713785 - 2.429260 N F 



-Nf 



3 
272 



192C(2)C(3) 



C(4) 



(4.8) 



Notice the appearance in N ~ l of three non-planar factors, namely C F (C F —Ca/2) and 
C f Ca{C f — Ca/2) and C F {C F —Ca/2)N f . These results are generalizing the observation of non- 
planarity of y)p ^"'-term of anomalous dimension function to three-loops and may be considered 



as the first non-obvious argument in favor of the correctness of definite results of Ref. [ [TI| ] 

Let us turn to the cot 
It has the following form 



Let us turn to the coefficient function of the valence contribution to the Gottfried sum rule P G . 



C%,) 



(l)N=l 



(4.9) 



Like in the case of the Adler sum rule, the first coefficient c[ N 1 is equal to zero (see e.g. 



Ref.[51]). The the second coefficients to the N-th Mellin moments of the NS combinations of 



structure functions F 2 of charged lepton-nucleon and neutrino-nucleon DIS are related to the Mellin 
moments from Eq.(2.11) as 



C, 



(l)N 



c 



(v)N 



3 ( l dx[C^ 
Jo 

i n 



dx\C {2) > 







x ) +C ^-\x)]x N - 1 
{x)-C^-\x)]^- x 



(4.10) 
(4.11) 
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In the case of N=l the explicit numerical integration of Eq.(llO) gave the following result [53] 



C 



(/) N=l 



3.695Q-1.847C A C F 



(4.12) 



where C^'^(x) in Eq.(fU0|) were taken from Ref.[^7] with the choice of MS renormalization 
scale jU 2 =<2 2 - Note, that non fulfillment of the second requirement, mentioned in the Introduction, 
namely the non-proper rounding off errors of numerical integrations of Ref.[53] did not allow 
to realize that Eq. ( 4.12 ) has the same non-planar structure as the expression for y[ 
fact was demonstrated in Ref. [Hj, where the following analytical result for C 



(l) N=l 



This 



(/) N=\ 



c 



(Z)AT=1 



(C 2 F -C A C F /2) 



141 21 45 

— + T C(2)-^?(3) + 12C(4) 



with the help of the Mellin representation 



C 



(l) N=l 



2x3 [ dx[C^-\x)\ 
Jo 



was obtained 
(4-13) 

(4.14) 



which follows from Eqs.(|4.5D, ( J4.6D and the requirement of absence of QCD corrections to the 
Adler sum rule. 

Next, taking the large N c limit for Eq.(4.1), which presumes that the corrections to the coeffi- 
cients of QCD j3-function of Eq.(2.3) and of Eq.(2.4) take the following form 



11 



,N F . 



(4.15) 
(4.16) 



and a s /{nN c ) ~ \/{N 2 ln{Q 2 /K 2 )), we find that Gottfried and Adler sum rules are rfelated as 



Tk 



h l + 0( 



Nl 



(4.17) 



This relation is scheme-independent. Indeed, the transformation of the a 4 -corrections in Eq. (|4JJ) 
to another M5-like scheme, can be done with the help of the shift 



a s {Q 2 ) _a s {Q 2 ) +fi Ja s {Q 2 ) x2 



% 



% 



K 



(4.18) 



where A is the concrete ^.-independent number, which is defined by the logarithm from the ratio 
of regularization scales and \^ s _i ike - Thus, the general MS-like scheme expression for the 
coefficient _1 takes the following form 

,(/)AT=l 



r V,iy-± _ (l)N=l 

^2 MS-like 2 MS 



(4.19) 



where both C^^LJ and Yi~ l contain the same non-planar group weight Cf(C f — Ca/2) [|37|]. 
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Keeping in mind the intriguing relation between the valence part of the Gottfried sum rule 
and the Adler sum rule, one may make a prediction, that the difference between NS moments of 
F2 in charged-lepton and neutrino scattering will continue to satisfy the property of non-planar 
suppression [^J. For example the ratio 



takes the value R% 1 



R 



latN 



C 



(l)N 



6C 



(v)N 



(XX 



N-l 



R 



■N=2 



1, but for N = 2, one cane obtained from Ref . J24]] 
0.50593104 



(4.20) 



(4.21) 



5.418324 liV? - 4N C N F - 8.4480127 

which is negative and small in magnitude at large N c , and is continuing to decrease rapidly with 
increasing of N (see discussions in Ref.[|6|]). So, we expect that the ratio 



6Q 



v)N 



(l)N 



1 + 0( a5> 



(4.22) 



is valid not only for 2 loops, but in higher orders of perturbation theory as well [Q]. 

It should be stressed, that the results and conjectures of Ref.[||], which were discussed in 
this talk, already received several confirmations. Indeed, in Ref. [|54|] the analytical expressions of 
Eq.Q and Eq.([U3]) were reproduced with the help of calculations, which are based on different 
theoretical technique. One more argument in favor of careful study of the topological structure of 
analytical perturbative series for characteristics of DIS |§] came from the recent analytical calcu- 
lations of the NS moments of vN DIS both at the cq and cr s level [55], and in the confirmation of 
the special role of non-planar-type colour structures in higher corrections of the quantities consid- 
ered. Thus, two works, presented at this workshop, satisfy the fourth requirement of importance 
of analytical calculations mentioned, in the Introduction, namely the importance of careful study 
of topological structute of perturbative QCD series. 



5. Large A^-limit: where are the effects of 0(N®) corrections to the Gottfried sum? 

First of all it it worth to stress that non-planar perturbative QCD corrections produce the 
following numerical results [H] 



1+0.035521^^ -0.58382^ 



1+0.038363) — 



. 0.56479 



for N F = 3 



for N F = 4 



(5.1) 



(5.2) 



which are higher than the results of different analysis of NMC data [|47|], though recently updated 
[18]. This discrepancy can be resolved by introducing to F the isospin-breaking effects in the 



Dirac sea u(x) < d(x) as in Eq.(3.5). It is interesting to note that within the large A^-limit of the 
modified soliton model of Ref. [|56|], proposed in Ref. [57], in the planar approximation one can 
get non-perturbative contribution of over 0(N®) [58] 



I! 



dx[ u{x)-d{x) =0(N^) 



(5.3) 
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which is absent in the perturbative series for P G . Thus, the 0(N~) -effects, which do not manifest 
themselves in perturbative sector, may reflect the existence of flavour-asymmetry terms in the nu- 
cleon sea. In the letter case the values of Iq between 0.219 and 0.178 were obtained for a range of 
constituents quark mass between 350 and 420 MeV, in fair agreement with the announced NMC 
result f* p = 0.235 ± 0.026 at Q 2 = 4 GeV 2 [0]. 

To summarize, we would like to emphasize that the bridge might exist between the 0(N®) 
corrections, which are absent in the perturbative sector, but are manifesting themselves in the the 
non-perturbative regime. 



6. Conclusions 

In this talk we would like to put extra attention on the importance of analytical multiloop an- 
alytical perturbative QCD calculations of characteristics of DIS. Indeed, the results of definite cal- 
culations reveal the existence of special topological structure in some physical quantities, related to 
phenomenology. Among them is the perturbative series for the Gottfried sum rule, which up to or- 
der a 2 corrections is defined by non-planar 0(1/ N 2 ) suppressed corrections only. It is interesting 
to understand, whether this typical non-planarity will continue in higher loops. The results of the 



talk of Ref. Q55p indicate that definite conjectures made in Ref. [Q] on continuation of non-planar 
suppression in the characteristics of vN DIS at the o? s level remain true. The possible physical 
explanation of this special structure of perturbative QCD series is the existence in the Dirac sea of 
isospin symmetry breaking effect, which is generating light quark flavour asymmetry. Another in- 
teresting physical feature of the similar origin is the perturbative generation of strange-antistrange 
asymmetry s(x) — s(x) in the nucleon sea ^ by NNLO corrections Np(Pq q — P s qq ) = a^P^g, where 
pjfe }' s term, calculated in Ref. [pT|] is proportional to d abc d a b c colour structure and result from "light- 
by-light" scattering type diagrams. Thus, one may expect, that the appearance of new colour 
structures in the perturbative series for DIS may be really connected to new physical conclusions, 
which should be studied more carefully. 
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